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ABSTRACT 

The Fermi Gamma-ray Space Telescope has revolutionized our knowledge of the gamma-ray pulsar 
population, leading to the discovery of almost 100 gamma-ray pulsars and dozens of gamma-ray 
millisecond pulsars (MSPs). Although the outer-gap model predicts different sites of emission for 
the radio and gamma-ray pulsars, until now all of the known gamma-ray MSPs have been visible in 
the radio. Here we report the discovery of a “radio-quiet” gamma-ray emitting MSP candidate by 
using Fermi , Chandra , Swift , and optical observations. The X-ray and gamma-ray properties of the 
source are consistent with known gamma-ray pulsars. We also found a 4.63-hr orbital period in optical 
and X-ray data. We suggest that the source is a black widow-like MSP with a ~ 0.1 M 0 late-type 
companion star. Based on the profile of the optical and X-ray light-curves, the companion star is 
believed to be heated by the pulsar while the X-ray emissions originate from pulsar magnetosphere 
and/or from intra-binary shock. No radio detection of the source has been reported yet and although 
no gamma-ray/radio pulsation has been found, we estimated that the spin period of the MSP is ~ 3 — 5 
ms based on the inferred gamma-ray luminosity. 

Subject headings: binaries: close - Gamma rays: stars - pulsars: general - stars: individual 

(1FGL J2339.7-0531, SDSS J233938. 74-053305. 2) - X-rays: stars 


1. INTRODUCTION 

The Fermi Large Area Telescope (LAT) has detected 
1873 y-ray point sources during its first 24 months of 
operation (Abdo et al. 2011). The majority of the extra- 
galactic sources have been identified as Active Galactic 
Nuclei (AGN; Abdo et al. 2010a), while many of the 
Galactic sources has been identified as 7-ray emitting 
pulsars (Abdo et al. 2010b). Gamma-ray pulsars are 
typically either young and energetic, such as the Crab 
and Vela pulsars, or very short period millisecond pulsars 
(MSPs). Throughout this paper we will refer to these two 
classes as the energetic and millisecond pulsars. In this 
second source catalog (2FGL), there are more than hun- 
dred identified energetic and millisecond pulsars. About 
half of them with gamma-ray pulsations were discovered 
previously by using radio pulsation search (e.g. Ransom 
et al. 2011 ; Keith et al. 2011; Cognard et al. 2011; Car- 
aveo 2010). In addition, although there were 35 gamma- 
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ray pulsars identified from the energetic class using a 
blind fourier search of the LAT data (see Pletsch et al. 
2011 and references therein), only four were subsequently 
identified as radio pulsars by folding the radio data with 
the pulsar ephemerides (Camilo et al. 2009; Saz Parkin- 
son et al. 2010; Pletsch et al. 2011). This indicates 
that there is a population of “radio-quiet” energetic pul- 
sars as seen from the Earth. However, all 21 gamma-ray 
emitting MSPs are “radio-loud”. Based on current ob- 
servations, we can classify the current population of 7- 
ray emitting pulsars into “radio-loud” energetic pulsars, 
“radio-quiet” pulsars and “radio-loud” MSPs. However, 
it is unclear if there is a class of “radio-quiet” 7-ray emit- 
ting MSPs. 

It is well known that the radio emissions of pulsars 
are associated with the activity of polar cap accelera- 
tor. On the other hand, the high-energy emissions from 
the pulsar magnetosphere have been studied with po- 
lar cap model (Ruderman & Sutherland 1975), slot gap 
model (Muslimov & Harding 2003) and outer gap model 
(Takata, Wang & Cheng 2010a). The polar cap model 
assumes the emission region is close to the stellar sur- 
face and above the polar cap, and therefore the model 
implies radio-loud 7-ray pulsars are much more common 
than radio-quiet 7-ray pulsars. The outer gap and slot 
gap models assume an acceleration region extending to 
the outer magnetosphere and indicates radio-quiet 7-ray 
pulsars are the major population of 7-ray pulsars. Based 
on the spectral shape in the GeV bands and the popula- 
tion of radio-quiet 7-ray pulsars found with Fermi (Abdo 
et al. 2010b), the outer gap or slot gap region is more fa- 
vorable as the origin of the 7-rays from energetic pulsars 
(Takata, Wang & Cheng 2010). 

For MSPs, “radio-quiet” millisecond 7-ray pulsars have 
not been identified so far. This may be due to following 
reasons: (1) the present sensitivity of blind frequency 
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search prevents a detection of millisecond pulsation in 
Fermi data; (2) most of known radio MSPs are in bi- 
naries that gamma-ray searches would be insensitive to, 
and (3) 7-ray emissions from MSPs always accompany 
radio emissions. On the other hand, Venter, Harding 
& Guillemot (2009) found that the pulse profiles of six 
MSPs detected by the Fermi can be fitted by the geome- 
tries with outer gap or the slot models. This implies that 
the emission regions of the 7-rays are different from the 
radio emission site for some of the MSPs and thus we ex- 
pect that there are “radio-quiet” 7-ray MSPs for which 
the radio beam is outside the line of sight. It is worth 
noting that “radio-quiet” here does not mean no intrinsic 
radio emission, just none beamed in our direction. 

In this Letter, we report a multiwavelength identifica- 
tion of a “radio-quiet” 7-ray MSP candidate. In a sepa- 
rate paper, Romani and Shaw (2011) arrived at similar 
conclusions via optical spectroscopy. 

2. UNIDENTIFIED FERMI SOURCE AS A “RADIO-QUIET” 
7-RAY EMITTING MILLISECOND PULSAR 

To identify suitable targets of “radio-quiet” 7-ray 
MSPs, we first selected candidates from the Fermi LAT 
first source catalog (1FGL; Abdo et al. 2010c) based on 
four criteria: 1) No known association at other wave- 
lengths; 2) source variability; 3) Galactic latitude, and 
4) gamma-ray spectral shape. 

We used the variability index in the 1FGL catalog to 
characterize source variability. Gamma-ray pulsars have 
always been found to be steady sources of gamma-ray 
emission (Abdo et al. 2010c, 2011). This property can 
therefore be used to help identify which of the uniden- 
tified sources are probably pulsars. The 1FGL catalog 
defines a variability index, for which a value greater than 
23.21 means that there is less than a 1% probability of 
being a steady source. We therefore selected sources with 
a variability index less than 23. 

MSPs are in general older than energetic 7-ray pulsars. 
Young objects like energetic 7-ray pulsars (r < 1 Myr) 
are likely located in the Galactic plane, while a frac- 
tion of MSPs should be at higher Galactic latitudes. We 
thus selected Fermi sources with high Galactic latitudes 
(|6| > 40°). The choice of this value is based on Monte- 
Carlo simulations for the Galactic population of MSPs, 
which show that no energetic 7-ray pulsar exists above 
|6| = 40° (Takata, Wang & Cheng 2011). 

Finally, we identified potential candidates from the 7- 
ray spectra. Although only a power-law spectrum is 
listed in the first-year catalog, it also has a curvature 
index to indicate how good of a power-law fit. For 7- 
ray pulsars, their 7-ray spectra is usually described by a 
power-law plus an exponential cutoff model (Abdo et al. 
2010b). According to the 1FGL catalog, the curvature 
index, for which a value greater than C = 11.34 indi- 
cates less than 1% chance that the power-law spectrum 
is a good fit. Hence we chose objects with a curvature 
index larger than 12. 

If a 1FGL source satisfied all four criteria, we short- 
listed as a potential candidate. From the short list, 
we further searched for X-ray imaging data ( Chandra , 
XMM-Newton and Swift ) from public archive and looked 
for X-ray sources within the error circles of Fermi as 
the first step for multiwavelength investigation. In this 
paper, we focus on one of our targets, 1FGLJ2339.7- 


0531 11 . 

3. MULTIWAVELENGTH IDENTIFICATION 

1FGL J2339. 7-0531 is one of the 205 bright gamma- 
ray sources detected with Fermi LAT during its first 
three months of operation (Abdo et al. 2009). It re- 
mains as an unidentified source in the second Fermi-L AT 
source (2FGL) catalog. In the 1FGL catalog (Abdo et 
al. 2010c), 1FGL J2339. 7-0531 has a variability index 
of 9.2 and a curvature index of 22.7. By comparing 
the 1FGL and 2FGL catalog (Abdo et al. 2011), the 
gamma-ray flux of 1FGL J2339. 7-0531 is constant and 
near 3 x 10 -11 ergs s -1 cm -2 . With better statistics, the 
second source catalog shows that a log parabola spec- 
trum can provide a better spectral fit comparing to a 
power-law model (Abdo et al. 2011). With its high 
Galactic latitude (—62°) and unidentified nature in all 
Fermi catalogs, 1FGLJ2339. 7-0531 is therefore one of 
the best targets to search for “radio-quiet” 7-ray MSPs. 

1FGLJ2339. 7-0531 was observed with Chandra and 
Swift/XRT for 21ks on 2009 October 13 (PI: Cheung), 
and 3.2ks on 2009 November 4, respectively. The Chan- 
dra imaging was operated with ACIS-I while the XRT 
was in the photon counting imaging mode. We repro- 
cessed with updated calibration files. HEASOFT version 
6 and CIAO version 4.2 were used for data reduction and 
analysis. Within the 95% error circle of Fermi (based 
on the refined position in the 2FGL catalog), there is 
only one relatively bright X-ray source in both obser- 
vations. The brightest X-ray source (CXOU J233938.7- 
053305) near the center of the error circle has an X-ray 
flux of 3 x 10 -13 ergs s -1 cm -2 (0.3-10 keV) based on 
a spectral fitting using an absorbed power-law with a 
best-fit photon index of 1.1 (based on the Chandra ob- 
servation). This indicates that the ratio between the X- 
ray and 7-ray (> 100 MeV) flux (F 7 ~ 3 x 10 -11 ergs 
s -1 cm -2 ) as measured by Fermi (Abdo et al. 2011) 
become Fx/F^ ~ 0.01, which is consistent with typi- 
cal observed values for 7-ray pulsars. It is worth not- 
ing that there are 9 much fainter Chandra sources in 
the error circle and their X-ray-to-gamma-ray flux ratios 
are less than 0.1%. Although we cannot totally rule out 
their association with the gamma-ray source, such a low 
flux ratio is not typical. We therefore identified CXOU 
J233938. 7-053305 as the potential X-ray counterpart to 
1FGLJ2339. 7-0531. 

Within the Chandra error circle (0.6 arcsec at 90% 
level) of CXOU J233938. 7-053305, there is a R ~ 19 
star from the USNO catalog (Monet et al. 2003). The 
same optical source is also seen in the SDSS Data Release 
8 images as SDSS J233938. 74-053305.2 with u' = 20.85, 
g' = 19.0, r' = 18.61, i! = 18.25, and z' = 18.23, as 
well as in the GALEX images with NUV=22.88 (177 
283 nm). Swift/UVOT observations taken simultane- 
ously with the X-ray observations detected the source 
with U = 19.58, B = 19.53, and V = 18.88. The X-ray- 
to-optical flux ratio (390) is much too high for a fore- 
ground star and an AGN (e.g. Green et al. 2004, Laird 
et al. 2009). Moreover, the spectral energy distribution 
from UV to optical indicates that it is unlikely to be an 

11 We will use the 1FGL naming conversion throughout the pa- 
per, and the source is called 0FGL J2339. 8-0530 and 2FGL J2339.6- 
0532 in the bright source catalog and 2FGL catalog, respectively. 
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Figure 1. Representative optical (white light) light-curve of the 
optical counterpart to 1FGL J2339. 7-0531, as observed by the 
Lulin lm telescope in Taiwan. The differential magnitudes are 
derived by comparison with several comparison stars in the field. 
Periodicity on a timescale of 4-5 hours is clearly seen. Also shown 
is the light-curve of a comparison star (with an average error of 
0.01-0.02 magnitudes). The timing resolution is 2 min. 

AGN (e.g. Richards et al. 2002). Instead, it looks like a 
F- or G-type star. In a separate paper (Romani & Shaw 
2011), they show via optical spectroscopy that the stel- 
lar object is consistent with a late-type star. We also 
had a long-term optical monitoring program with the 
MITSuME 50cm telescope (Kotani et al. 2010) located 
at Akeno Observatory in Japan from 2010 September to 
2010 November. MITSuME equips with a tricolor cam- 
era that can perform simultaneous imaging in the g ' , R , 
and / bands. It is clear from Table 1 that all three bands 
show noticeable variability. Based on the above obser- 
vations, we believe that the gamma-ray/X-ray source is 
likely a binary system with optical emission from the 
irradiation of the companion star, while the central com- 
pact object is responsible to the X-ray and gamma-ray 
emission. 

Because 1FGL J2339.7-0531/CXOU J233938.7 

053305 may be a new class of interesting object, we 
carried out an intensive optical monitoring campaign 
using the lm telescope at the Lulin Observatory in 
Taiwan and the 0.8m Tenagra Telescope in Arizona 
(see Table 1 for the observation log). All images are 
flat-field and bias corrected and we performed relative 
photometry by comparing with several comparison stars 
in the field. In all 5 nights, the light-curves clearly show 
variability on a timescale of 4-5 hours (see Figure 1 
& 2). We then performed a timing analysis by using 
the Lomb-Scargle periodogram on the combined optical 
data and found a period of 4.6342(9) hours. Figure 
2 shows the folded light-curves of all the optical data 
(upper panel) and Chandra data (lower panel) at a 
period of 4.6342 hours. The phase zero is chosen at 
MJD 55500. In addition, we folded the optical colors 
( g ' — I) obtained with MITSuME and X-ray hardness 
ratio (1.5-8keV/0.3-1.5keV) at the same period (Figure 
2). 

4. DISCUSSION 

Using multi-wavelength data, we found an X- 
ray/optical counterpart to 1FGL J2339. 7-0531. The X- 
ray-to-gamma-ray flux ratio is consistent with a gamma- 
ray pulsar while the gamma-ray spectrum deviates from 
a power-law which is commonly seen in MSPs. The X- 
ray spectrum is very hard with a power-law photon index 


of 1.1. Such a hard spectrum is not typical for X-ray bi- 
nary, AGN, and foreground star. However, some MSPs 
show similar spectral characteristics (e.g., Archibald et 
al. 2010; Tam et al. 2010) and such a hard spectrum 
could be from a pulsar wind nebula. We propose that 
1FGL J2339. 7-0531 is a gamma-ray emitting MSP in a 
binary system and the optical emission is from a late- 
type companion star with contribution from the heating 
on the stellar surface. 

The optical light-curve shows clear modulation at a pe- 
riod of 4.6342 hours (Figure 2). It is beyond any doubt 
that the periodicity is associated with the orbital period 
of the system. Given the large variation (about 3-4 mag- 
nitudes) in the optical light-curve, the inclination must 
be high. Better photometric data in the near infrared 
and light-curve modeling will be able to constrain the 
inclination in the future. Assuming the companion has 
filled its the Roche lobe, the mass limit of the companion 
star can be estimated from the orbital period. For a 4.6- 
hr orbital period, the mass for a normal main-sequence 
star is < 0.5 M 0 (Wilson et al. 1999). 

The large optical variability and hints of color vari- 
ation (Fig. 2) suggest that the companion is being 
heated by the pulsar. The irradiation of pulsar 7-ray 
emissions onto the companion star can produce the or- 
bital modulation of the optical emission from the bi- 
nary system (e.g. Takata, Cheng & Taam 2010, 2011). 
In this model, the optical maximum occurs at the in- 
ferior conjunction, where the pulsar is located between 
the companion star and the Earth. Furthermore, dur- 
ing the optical maximum, the color of the star tends 
to be bluer because of the heating effect. There is 
some indication from the MITSuME data that the color 
of the star has become bluer during the optical maxi- 
mum (Fig. 2). The maximum luminosity of the optical 
emission can be estimated as SL opt ~ {-kO 2 / 8Vl)L 1 ~ 
lO 31 (0/O.lrad) 2 (JU/3rad) _1 (L 7 /lO 33 ergs s _1 ) ergs s _1 , 
where 6 is the angular size of the companion star mea- 
sured from the pulsar and 8fl is the solid angle of the 
7-ray beam. The observed maximum R-band magnitude 
of 1FGL J2339. 7-0531 is about 18 (see Table 1) which is 
consistent with the above estimation. 

It has been suggested that X-ray emissions from black 
widow systems originate from pulsar magnetosphere 
and/or from intra-binary shock due to the interaction 
between the pulsar wind and the injected material from 
the low mass star (e.g., Kulkarni et al. 1992, Stappers 
et al. 2003). If the X-rays are produced in the magneto- 
sphere, the orbital modulation will be produced due to 
the physical eclipse of the pulsar by the companion star. 
In this case, the minima of the X-ray and optical light- 
curves should be at the same phase, and we indeed ob- 
served that (see Fig. 2). Alternatively, if the X-rays orig- 
inate from the intra-binary shock region, an effect of the 
Doppler boosting relating with the post-shock flow can 
produce the orbital modulation (Arons & Tavani 1993). 
It has been suggested for the pulsar and massive stellar 
system (e.g. PSR B1259-63/LS 2883 system, Bogovalov 
et al. 2008), the post-shock flow can be accelerated into 
relativistic regime because of an explosion in the down- 
stream region (see also Tam et al. 2011). The Doppler 
boosting increases the intensity as I v oc D 3+a I„, where 
is the intensity in the comoving frame, D is Doppler 
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Table 1 

Optical Observation log of 1FGL J2339. 7-0531 


Date and Time (UT) 

Telescope 

Filter 


Exposure 0, 

2010-10-25 

14:44-15:57 

Lulin lm 

V 


2 min 

2010-10-26 

13:33-15:59 

Lulin lm 

V 


2 min 

2010-10-31 

10:23-15:53 

Lulin lm 

White 


2 min 

2010-11-01 

11:33-16:27 

Lulin lm 

White 


2 min 

2010-11-11 

01:45-06:20 

Tenagra 0.8m 

White 


5 min 

2007-09-13 

22:32:08 

GALEX 

NUV=22.88 


1642 sec 

2008-11-20 

05:25-05:29 

SDSS u' = 20.85, g' = 19.0, r' = 18.61, 

i' = 18.25, z' = 18.23 

54 sec 

2009-11-04 

07:56-16:07 

UVOT 

U = 19.58, B = 19.53, V = 18.88 

1144, 838, 1144 sec 

2010-09-13 

14:52:26 

MITSuME 0.5m 

g' = 18.51, R = 18.04, 

/ = 17.53 

60 sec 6 

2010-09-17 

14:46:37 

MITSuME 0.5m 

g' = 18.54, R = 17.98, 

/ = 17.63 


2010-09-20 

13:17:30 

MITSuME 0.5m 

g' = 18.15, R = 17.33, 

I = 17.38 


2010-10-04 

15:55:20 

MITSuME 0.5m 

g' = 18.51, R = 17.93, 

/ = 17.53 


2010-10-10 

14:28:06 

MITSeME 0.5m 

g' = 19.12, R = 18.53, 

I = 17.94 


2010-11-09 

13:21:22 

MITSeME 0.5m 

g' = 18.43, R = 17.82, 

I = 17.21 


2010-11-10 

13:22:43 

MITSeME 0.5m 

g' = 18.98, R = 18.12, 

I = 17.51 


2010-11-11 

13:19:40 

MITSeME 0.5m 

g' > 20.19, R = 19.44, 

I = 19.04 


2010-11-12 

11:22:28 

MITSeME 0.5m 

g' = 18.43, R = 17.70, 

/ = 17.44 


2010-11-14 

12:39:55 

MITSeME 0.5m 

g' = 19.00, R = 18.16, 

I = 17.84 


2010-11-15 

12:40:11 

MITSeME 0.5m 

g' = 18.07, R = 17.91, 

I = 17.51 


2010-11-18 

12:53:05 

MITSeME 0.5m 

g' > 20.02, R = 21.77, 

I = 19.16 


2010-11-21 

12:46:57 

MITSeME 0.5m 

g' = 18.70, R = 17.61, 

I = 17.49 



a Exposure time for each frame. 

b Some frames were combined to increase the signal-to-noise ratio. 


factor and a is the spectral index in the co-moving frame. 
An orbital modulation by a factor of five may suggest the 
Doppler factor D ~ 1.6 with a ~ 0.5. 

The X-ray hardness ratio roughly correlates with the 
X-ray light-curve (Figure 2). It may suggest that the pul- 
sar (which has a very hard X-ray spectrum) is physically 
eclipsed by the companion and that the soft emission 
is an extra source of light (e.g. a wind shock). There 
is some support for this in the assymetry of the optical 
light curve. Much better X-ray and optical photometric 
data are required for serious modeling. 

The source has a very low absorption (3 x 10 20 cm -2 ) 
along the line of sight. In comparison with Geminga with 
Nh = 1.5 x 10 20 cm -2 and a distance of 150-250 pc, 
and SAX J1808. 4-3658 with Nh = 4 x 10 20 cm -2 and a 
distance of 1.3 kpc. The distance of 1FGL J2339. 7-0531 
is therefore likely to be between 300 and 1000 pc. If the 
companion star is a late-type main sequence (e.g., a M5 
dwarf), and has V ~ 21 during optical minimum when 
the heating effect is at minimum, the implied distance is 
about 700 pc. 

If we assume a distance of 700 pc, then the gamma-ray 
power is i 7 ~ 10 33 ergs s _1 . Most of MSPs have an 
efficiency about 10 % which gives spin-down luminosity 
L s d ~ 10 34 ergs s . According to Takata, Cheng & 
Taam (2011), the spin period of the MSP can be derived 
from the spin-down power. Using the estimated spin- 
down power, the spin period will be ~ 4 ms. We note 
that the distance to the source is uncertain; for instance, 
Romani and Shaw (2011) estimated a distance of 1.1±0.3 
kpc. If we assume a distance error of a factor of 2 (i.e. , 
350-1400 pc), the spin period will be 3-5 ms. 

The bright Fermi source 1FGL J2339. 7-0531 is an 
intriguing object that is potentially the first “radio- 
quiet” 7 -ray emitting MSP in a binary system. Despite 


1FGL J2339. 7-0531 being one of the brightest unidenti- 
fied pulsar-like Fermi sources located at high Galactic 
latitude, no radio detection has been reported. Based 
on a recent pulsation search with the Green Bank Tele- 
scope, 1FGLJ2339. 7-0531 was not detected (Ransom et 
al. 2011). It may indicate that 1FGL J2339. 7-0531 is 
“radio-quiet” , or the radio emission must be very weak. 

The system resembles black widow MSP for which the 
MSP has a very low-mass (<< O.IMq) companion in 
an orbit less than a day. The optical light-curves of 
black widow systems usually shows large orbital varia- 
tion because the irradiation produces strong heating on 
the companion facing the pulsar. It is therefore sugges- 
tive that the companion (usually a white dwarf) is being 
evaporated by the high-energy radiation from the pul- 
sar. On the other hand, there is a class of black widow 
MSPs with high mass (> O.IMq) non-degenerate com- 
panions. One classic example is PSR J1023+0038 that 
has a 0.2Mq companion (Archibald et al. 2009). The 
X-ray /gamma-ray properties of 1FGL J2339. 7-0531 are 
similar to that of PSR J1023+0038 but no radio emission 
from the former has been detected yet. 

Hence, we suggest that 1FGL J2339. 7-0531 is a new 
class of black widow-like MSP system, with no or 
very faint radio emission along the line-of-sight. The 
high-mass non-degenerate companion indicates that 
1FGL J2339. 7-0531 is in the late stages of recycling. The 
absence of radio emission is the result of different emis- 
sion regions for radio and gamma-rays. For instance, the 
radio emission is from the polar cap with the radio beam 
out of the line-of-sight. On the other hand, the gamma- 
ray emission is from the outer magnetosphere which is 
predicted by the outer gap model. 

An intensive search for radio emission and pulsation 
from 1FGL J2339. 7-0531 in the future can confirm the 
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Figure 2. Folded light-curve of optical and Chandra observations of 1FGL J2339.7— 0531 with a best-fit period of 4.6342 hr. Optical colors 
(g' — I) obtained with MITSuME are also plotted. The phase zero is defined as 2010 October 31 (MJD 55500). Also plotted with the 
X-ray light-curve is the X-ray hardness ratio (1.5—8 keV/0.3-1.5 keV) with triangles. It is evident that both optical and X-ray light-curve 
show similar modulation. The X-ray hardness ratio exhibits some variability when the X-ray light-curve is at its minimum. Note that the 
Chandra data only cover about one orbital period. 


above scenario. At the same time, a search for gamma- 
ray/X-ray pulsation from the MSP will reveal the true 
nature of the system. 
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